Cosmological Solution to Einstein- Vlasov System 

Tigran Aivazian ^] 
Abstract 

Einstein- Vlasov system is solved for a homogeneous isotropic spacetime with positive 
curvature (R x S 3 topology). For the Universe consisting of massless particles the equation 
for R(t) is solved analytically. 



1. General Equations 

We start with well-known system of self-consistent field equations of Einstein- 
Vlasov which, in geometrical form can be written as: 

Ein(g) — J fp®pd(x (1.1) 

Xf = (1.2) 

where X is a vector field which generates phase flow. We use geometrical 
system of units where c = 1,G = 1/87T. In coordinate form the system (1.1) 
can be rewritten as: 

- l/2R 9flu = -7=1 f(x,p)p^ Pu 6(p 2 + m 2 )d 4 p (1.3) 
v Q 

9 p ^ = 2^%: (L4) 

where f(x,p) is a distribution function of matter consisting of particles with 
mass m > 0. It is assumed to be a non-negative function with compact support 
on a mass shell in cotangent bundle T*M: 

S m (M) = {(x,p) e T*M I g(p,p) +m 2 = 0}, 

/ : S m (M) -> R + , 
S m (M) = Si(M)US^(M) } 

F± = f I si (1.5) 

Note, that our expression for T^ix) differs from the one used in [1] only slightly, 
which allows to avoid Christoffel symbols in the equations for characteristics. 
This corresponds to working in (x^,p u ) coordinates as opposed to (x M ,p v ). 

It is convenient to use ADM 3+1 [2] representation to solve Cauchy problem. 
In Appendix A we provide the 3+1 representation of Einstein- Vlasov system 
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for a generic metric. In synchronous frame of reference, metric tensor has the 
following form: 

^ = G L) (L6) 

The dynamical equations are as follows: 

d tlik = -2K ik (1.7) 

d t K ik = R lk - 2K U K\ + KK ik + ]p lk ^ lm T lm - T 00 ) - T lk (1.8) 
Constraint equations are: 

R + K 2 -K ik K ik = 2T 00 (1.9) 

D j {K\-5\K) = -T 0i (1.10) 

where i? = Y k Kik, Di is a covariant derivative in 7^ metric. It is not difficult 
to see that kinetic equation for f(x,p) will split into a couple of equations 
for F + (x,p,t) and F_(x,p,t) which are constructed by reducing f(x,p) on 
"upper" and "lower" halves of mass shell. The reason for being able to reduce 
an equation on S m (M) into a pair of equations for each connected half of it is 
due to the fact of S m (M) is dynamically invariant with respect to the phase 
flow generated by the vector field X. Here, we changed the meaning of what is 
implied by (x,p) as it is usually trivial from context to guess to which space 
it belongs to. Explicit form of kinetic equations becomes: 

^ ± .__( yV _-_^ M _) f± (L11) 

or, in terms of 3-Hamiltonian: 

p%d t F± = {H,F±} (1.12) 

H(x,p,t) = ^ ik (x,t) PiPk (1.13) 

Presence of p in these equations points to the fact that equations of character- 
istics for reduced kinetic euqations do not necessarily constitute a Hamiltonian 
system. The metric of the spacetime can be written in Robertson- Walker form: 

ds 2 = -dt 2 + R 2 {t)[d X 2 + sin 2 X (de 2 + sin V# 2 )] (1-14) 
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Then it is not hard to calculate the components of energy-momentum tensor 
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k = 



R 3 sin 2 x sin 9 ' 



Too = k 



r {F + + F_) PiPk < 
J ^ m 2 + hM 




T 0l = kJ (F+ - F_) Pl d 3 p 



(1.15) 



where X(p,p) = pi + p\j sin 2 x + pl/(sin 2 x sin 2 9). By the following transfor- 
mation A is brought to canonical form: 



In k-coordinates it is not hard to see the set of solutions for kinetic equations: 



where k = kf + k 2 + k 2 . The first constraint equation takes the following form: 



where we denoted F(k) = F + (k) + F_{k). It is interesting to notice that 
having set a restriction on the topology of the spacetime we cannot distinguish 
between particles and anti-particles, keeping in mind Feynman's Principle of 
Reinterpretation and interpreting values of f(x,p) in the area po < ("lower" 
half of mass shell) as the distribution function of anti-matter. On one hand, 
having F_ in energy-momentum tensor may violate dominant energy condition 
but, on the other hand, it appears to be more consistent and complete. Nature 
is symmetric and beautiful and so must be the equations that describe it 
however naive this may sound. One can derive second order equation for R(t) 



h = pi, 
k 2 =p 2 / sinx, 
h = P3/(sinxsin0) 



(1.16) 



F± = F ± (k) 



(1.17) 




(1.18) 



from (1.7) and (1.8): 




(1.19) 
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It can be easily shown that the constraint equation (1.10) is satisfied. Also, 
it is obvious (in this particular case) that if R(t) satisfies (1.18) for arbitrary 
t then it satisfies (1.19) as well which allows one to consider (1.18) as a main 
dynamical equation. The "maximal" set of solutions of (1.19) is, however, 
larger than that of (1.18) but we shall not be concerned with losing some of 
the solutions for now. 



2. Analytical Solution 

In this section we consider a special case of a Universe consisting of massless 
m=0 particles. Using dimensionless variables (x,r) we have: 

R = R(0)x, 
t = R(0)t, 
dx 

e = ii?(0) 2 P (0) (2.1) 
Equation (1.17) then becomes: 

x 2 = -i + 4, 

X A 

x(0) = 1 (2.2) 
If e > 1 Cauchy problem (2.2) has a real solution: 

x(t) =v / e - ( r -v^T) 2 (2.3) 
Going back to (t,R) variables we have: 



R(t)=Ri0 ^mm. (2 . 4) 

Using the above expression one can evaluate the age of the Universe (T^) and 
current value of the Hubble constant: 

R(0) ~ 1.25 • 10 28 cm, 
p(0) ~ 1.1 • 10~ 57 cm" 2 , 
e ~ 1.459, 
T bb ~ -7.09 • lOVs, 
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T tot ~ 31.7 • 10 B yrs, 
H = §^ ~ 5.304 • KT'W 1 (2.5) 

As follows from these calculations even under assumption m=0 the results are 
reasonable ie they lie within the " standard" bounds predicted by hydrodynam- 
ical Friedmann-Robertson- Walker cosmological model. More complicated case 
of m > can serve as a topic for next paper on the subject. 

Acknowledgements. 

This work was done several years ago (1992) when author was an undergrad- 
uate student of Yerevan State University (Armenia) and it remained in the 
form of "half-forgotten manuscript" until recently I read an article by Rendall 
[1] which woke up my interest to the Einstein- Vlasov system and incentive to 
"dig out old papers" . 

References. 

[1] A. D. Rendall.: An introduction to the Einstein- Vlasov system. |gr-qc/9604001 



[2] Arnowitt R., Deser S., Misner C.W., in Gravitation: Introduction to Cur- 
rent Research, ed. L. Witten, New York, Wiley, 1962. 



5 



3. Appendix A: ADM Representation of Einstein- Vlasov system 



Metric tensor: 

-a 2 + (3i(3 l (3. 



*» = < ft ;"j (31) 



Energy-momentum tensor: 
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T^(x) = -— / f(x,p)p,p u 5(p 2 + m 2 )rf 4 p (3.2) 
v fi 1 

Measure of integration in the phase space: 

dfi= jL=5(p 2 + m 2 )d 4 p (3.3) 

Dynamical equations (c = 1, G = 1/Stt): 

d tlij = -2aK l3 + D i (3 j + Djfa (3.4) 

d t Kij = —DiDjOi + a(R i3 - 2K u K l 3 + KK i3 - S i3 + 

\lrAl kl Tki - p m \) + + KuDj/3 1 + KtjDi/3 1 (3.5) 

d t F ± = L ± F ± (3.6) 

Constraint equations: 

R + K 2 - K lk K ik = 2p m (3.7) 

Dj(Kj - 5i j K) = Si (3.8) 

The expression for Liouville operators can be obtained from the equations for 
geodesies: 

2 t^* + (39) 

We use the following standard notations: 

F± = f{x,p) \ po=p ±, 
n^-a, 0,0,0), 
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ra"(l/a,-/37c0, 

ym - L fj,u 11 ' ,u 5 
Si T^yTl^P j, 

C 'f pv 

^ik - 1 - fiv 1 i 1 ki 

P\ = b\ + n»n v , 
K = Y k K ik , 
R = 7* fe .Rife, 

p l l I l m m I 

^ik ; ik,l 1 il,k ~r 7 ik] Im ! ill km,) 

= —^-{ljm,k + lkm,j ~ ljk,m) (3.10) 
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